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Abstract 

We propose a method for computing any Gelfand-Dickey r function living in Segal- 
Wilson Grassmannian as the asymptotics of block Toeplitz determinant associated to a 
certain class of symbols W(t;z). Also truncated block Toeplitz determinants associated 
to the same symbols are shown to be r function for rational reductions of KP. 
Connection with Ricmann-Hilbcrt problems is investigated both from the point of view 
of integrable systems and block Toeplitz operator theory. 
Examples of applications to algebro-geometric solutions are given. 



Introduction 



This paper deals with the applications of block Toeplitz determinants and their asymptotics 
to the study of integrable hierarchies. Asymptotics of block Toeplitz determinants and their 
applications to physics is a developing field of research; in recent years it has been shown 
how to compute some physically relevant quantities (e.g. correlation functions) studying 
asymptotics of some block Toeplitz determinants (see [25j,[26j,[27J). In particular in [25J and 
|26] authors, for the first time, showed effective computations for the case of block Toeplitz 
determinants with symbols that do not have half truncated Fourier series. This is of particular 
interest for us as, with our approach, we will be able to do the same for certain block Toeplitz 
determinants associated to algebro-geometric solutions of Gelfand Dickey hierarchies. Let 
us mention some theoretical results about (block) Toeplitz determinants we will use in this 
paper. Given a function 7(2) on the circle we denote Tn{j) the Toeplitz matrix with symbol 
7 given by 



(0) 



7 



(-AT) 



\ 



7 



(-AT+1) 



7 (0) J 

where are the Fourier coefficients 7(21) = Y^k 1^ z k ■ 

We use the term block Toeplitz for the case of matrix- valued symbol 7(2). In that case the 
entries j^~" i > of the above matrix are n x n matrices themselves. We denote 

D N (j) : = detT/v(7) 

and we use the notation T(7) for the N x N matrix obtained letting N go to infinity. The 
main goal of the theory of Toeplitz determinants is to compute Dn{j) as N goes to infinity 
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and find expressions for Djv(t) as well as for its limit in terms of Fredholm determinants. 
First result is due to Szego that in 1952 gave a formula for asymptotics of Dfffa) in the scalar 
case [5]. This result has been generalized by H. Widom in the 70's ([6], [7] and [8]) for the 
matrix case; namely he proved that under suitable analytical assumptions it exists the limit 

DM := ^oo^ =det(T(7)r( ^ 1)) 

where G(7) is a normalizing constant and the operator T(7)T(7 _1 ) is such that its determi- 
nant is well defined as a Fredholm determinant (see section 2 for the precise statement). Once 
the asymptotics had been computed the next quite natural question was to find an expressions 
relating directly D^^j), and not just its asymptotics, to certain Fredholm determinants. The 
problem was solved many years later by Borodin and Okounkov in [9] for the scalar case and 
generalized, in the same year, for matrix case by E. Basor and H. Widom in [10]. For matrix 
valued case Borodin-Okounkov formula reads 

D N ( 1 ) = D OQ ( 7 ) det(I - iCy.jv) 

(here we assume G{^) = 1). The operator (/ — if-^jv) can be written explicitely in coordi- 
nates knowing certain Riemann-Hilbert factorizations of 7. Its Fredholm determinant is well 
defined (see section 2 for details). Now many proofs of Borodin-Okounkov formula are known 
(for instance [TT] contains another proof of the same formula, see also the earlier paper [12j). 



In this paper we apply block- Toeplitz determinants to the computation of r function of 
an (almost) arbitrary solution of Gelfand-Dickey hierarchy 



(L differential operator of order n, j 7^ nk). 
More precisely to a given point 

W = W{z)H { + ] 



in the big cell of Segal- Wilson vector- valued Grassmannian we associate anxn matrix- valued 
symbol W(i; z) obtained deforming W(z) (see formula f 1 1 f) ) . In this way we define a sequence 
of ./V-truncated block Toeplitz determinants {Tw,N(t)}N>o which are shown to be solutions 
of certain rational reductions of KP; this is our 

First result: Every symbol W(i; z) defines through its truncated determinants a sequence 
{T\V,N(t)}N>o of solutions for KP such that 

r W ,N(t) E cKPi in7V n cKP n , n ViV > 0. 



Here we used the notation from [20] : given a r function for KP with corresponding Lax 
pseudodifferential operator C we say that r 6 cKP mjra iff C m can be written as the ratio of 
two differential operators of order m + n and n respectively. This sequence admits a stable 
limit which is shown to be equal to the Gelfand-Dickey r function r^y(t) associated to W; this 
quantity can be computed using Szego- Widom's theorem. This will give us the remarkable 
identity 

Twit) = det V w(i . z) (1) 
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where V W (f-z) * s ^ e Fredholm operator appearing in Szego-Widom's theorem (here we put t 
instead of t to remember that, when working with W € Gr^ n \ times t n j multiple of n must be 
set to 0). Next step is the study of Riemann-Hilbert (also called Wiener-Hopf) factorization 
of symbol W(i; z) given by 

W(t;z) =T-{t;z)T+(t;z) (2) 
with T_ and T+ analytical in z outside and inside S 1 respectively and normalized as 

T_(oo) = I. 

Here we assume that the symbol can be extended to an analytic function in a neighborhood 
of S . Using Plemelj's work [T3] we show that T—(t;z) must satify the integral equation 

rw&.) T -& z ) = I ( 3 ) 

and we write a solution of (|2J) in terms of wave function ipw{t) z ) corresponding to W. 
In this way we arrive to our second result: 

Second result: Take W £ Gr^ n ' in the big cell and its corresponding r function Tyy(i). 
rw(t) is equal to the Fredholm determinant of the homogeneous integral equation associated 
to (Ep which is related to Riemann-Hilbert problem (Q|). The solution of this Riemann-Hilbert 
problem is unique for every value of parameters t that makes Tw{i) / and can be computed 
by means of related wave function ipw{t; z). 

At the end of the paper we consider a particular class of symbols W(t; z) corresponding 
to algebro- geometric solutions of Gelfand-Dickey hierarchies. We formulate an alternative 
Riemann-Hilbert problem equivalent to © and explain how to solve it using ^-functions. In 
this way we give concrete formulas for a wide class of symbols that do not have half truncated 
Fourier series. We think this is quite remarkable since concrete results for non half truncated 
symbols were available, till now, just for the concrete cases presented in [25] and |26j. 
The paper is organized as follows: 

• First section states some results about Segal- Wilson Grassmannian and related loop 
groups we will need in the sequel; proofs can be found in pQ and [2]. 

• Second section states Szego-Widom's theorem and related results obtained by Widom 
in [6], [7] and [8] and the Borodin-Okounkov formula for block Toeplitz determinant |10j . 

• In the third section we introduce and study the sequence of truncated determinants 
{T~w,N(t)}N>o an d its stable limit Twit)- We want to remark that the property of 
stability was stated for the first time in [pS] (see also [TB] ) and our sequence is actually a 
subsequence of the stabilizing chain studied in [T7]; nevertheless, to our best knowledge, 
this is the first time that block Toeplitz determinants enter the game and also the 
observation that tw,n £ cKP n .„ seems to be something new. 

• Fourth section is devoted to establishing the connection between integral equations for- 
mulated by Plemelj in [13] and Fredholm operator appearing in Szego-Widom's theorem. 

• In the fifth section we show how to write Riemann-Hilbert factorization of W(t; z) in 
terms of wave function ip\y{t;z). Of course relation between Gelfand-Dickey hierarchy 
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and factorization problem is something known; our exposition here is closely related 
to [2]. Moreover, knowing Riemann-Hilbert factorization of W(t;z), we can apply 
Borodin-Okounkov formula to give an expression of any Tw,N{i) as Predholm determi- 
nant and a recursion relation to go from Tw,N(t) to Tw,N+i\t)- 

• Last section gives explicit formulas for symbols and r functions associated to algebro- 
geometric rank one solutions of Gelfand-Dickey hierarchies. Also we formulate an alter- 
native Riemann-Hilbert problem equivalent to ([2]) in analogy with what has been done 
in [25] and [26]. We explain how to solve it using ^-functions. 
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1 Segal-Wilson Grassmannian and related loop groups 

Here we recall some definitions and results from [1] and [2] that will be useful in the sequel. 

Definition 1.1. Let := L 2 {S 1 ,C n ) be the space of complex vector-valued square-integrable 
functions. We choose a orthonormal basis given by 

{e aM := (0,...,z fc ,...,0) T :a = l...n,keZ} 

and the polarization 

H (n) = R (n) R (n) 

where and are the closed subspaces spanned by elements {e a ^} with k > and 

k < respectively. 

In the sequel in order to avoid cumbersome notations we will write H instead of . 

Definition 1.2 Q2J). The Grassmannian Gv(H^) modeled on H^ n > consists of the subset of 
closed subspaces W C ij( n ) such that: 

(n) 

• the orthogonal projection pr + : W — ► H\ is a Fredholm operator. 

(n) 

• the orthogonal projection pr_ : W — > H_ is a Hilbert- Schmidt operator. 
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Moreover we will denote Gr^ the subset of Gr(H^) given by subspaces W such that zW C 
W. 



It's well known pQ that through Segal- Wilson theory we can associate a solution of 



n 



tli 



Gelfand-Dickey hierarchy to every element of Gr^ n ^ ; this is the reason why we are interested 
in them. 

Lemma 1.3 ([I]). The map 

E : — ► H 

(/„(*), f n ^(z)) T /(*) := f (z n ) + ... + ^-Vn-i(^) 
is an isometry. Its inverse is given by 

aw = \ E c~ fe /(o 

where the sum runs over the n th roots of z. 

Proposition 1.4. Under the isometry H we can identify Gn n ) with the subset 

{W e Gt(H) : z n W C W} 

It is obvious that loop groups act on Hilbert spaces defined above by multiplications. 
We want to define a certain loop group LiGlfn, C) with good analytical properties acting 

2 

transitively on Gr^; in such a way we can obtain any W £ Gr^ n ^ just acting on the reference 

(n) 

point H\ with this group. Good analitycal properties will be necessary as we want to 
construct symbols of some Toeplitz operators out of elements of this group and then apply 
Widom's results (see below). Given a matrix g we denote with \\g\\ its Hilbert-Schmidt norm 

n 

II Il2 n 1 1 2 

Definition 1.5. Given a measurable matrix-valued loop 7 we define two norms 1 1 7 1 1 00 and 

II7IL 1 as 
11 uz, 2 

1 

\2\ 2 



:= ess sup || 7 (*)|| \\j\\ 2j i := ^ ( ||7 (fc) || 



where we have Fourier expansion 

00 



(*)** 

k=— 00 



Definition 1.6. LiGl(n,C) is defined as the loop group of invertible measurable loops 7 such 

2 

that 

IMloo + H7II2 I < °°- 

Proposition 1.7 ([2j). LiGl(n, C) acts transitively on Gr^ and the isotropy group of 
is the group of constant loops Gl(n, C). 
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Proof can be found in [2], here we just mention the principal steps necessary to arrive to 
this result. 

• We define a subgroup G\ res (H^) of invertible linear maps g : 

R (n) _^ R (n) acti 

on 

Gr(H^) (the restricted general linear group). 

• We prove that every element of Gl res (H^) commuting with multiplication by z must 
belong to LiGl(n,C). 

2 

• We take an element W £ Gr^ ra ^ and a basis {w±, ...,w n } of the orthogonal complement 
of zW in W. 

• Out of this basis, putting vectors side by side, we construct W and easily check that 
W = W{z)H { ^ ) . 

• We verify that multiplication by W belongs to Gl res (H^); since it obviously commutes 
with multiplication by z we conclude that W(z) G LiGl(n,C). 



2 Szego-Widom theorem for block Toeplitz determinants. 

In his work ([6], [7] and [5]) H. Widom expressed the limit, for the size going to infinity, of 
certain block Toeplitz determinants as Fredholm determinants of an operator V acting on 
H^. Also he gave two different corollaries that allow us to compute this determinant in 
some particular cases. In this section we recall, without proofs, these results. Moreover we 
state Borodin-Okounkov formula as presented in [10J for matrix case. 

We begin with some notations; given a loop 7 G LiGlfn, C) we denote with T^Oy) the block 

2 

Toeplitz matrix given by 



(1 



(o) 



7 



(-7V) 



\ 



7 



(-AT+1) 



\7 W 



v(0) 



where we have the Fourier expansion 7(2) = ^ fc 7^'^z fc .We denote D^ipf) its determinant. 
We use the notation T(-f) for the N x N matrix obtained letting N go to infinity. 

Remark 2.1. It's easy to see that, in the base we have chosen above for H^ n \ T(j) is nothing 
but the matrix representation of 



pr + o 7 : — ► H\ 
Theorem 2.2 (Szego-Widom theorem, [5]). Suppose 7 E LiGl(n, C) and 



r(«) 



A arg det (7(e l °) 

O<0<2tt V V 
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Then it exists the limit 



where 



1 



2- 



G{i) = exp [— J log {det 7 (e ie ))de 

The proof of the theorem is contained in [8]; instead of rewriting it we simply consider 
the operator T(7)T(7 _1 ) and explain the meaning of "det" in this case. 

Lemma 2.3. Consider 71, 72 £ LiGl n (n,C); we have 



T(7i7 2 )-n7i)n7 2 ) = [E^ +/ " u_ " /,) 



fc>l 



Proof The (i, j)-entry of left hand side reads 



E 7?-*^ - E^r k H k ' J) = E 7r fc) 7^- J) = E7i' 

fc=— 00 fc=0 fc=oo fc=0 



In particular choosing 71 = 7 and 72 = 7 1 we obtain 

/-T(7)T(7- 1 ) = [^7^)( 7 -i)(-^fc) 



(i+k+i) (-k-j-i) 

72 



□ 



k>l 



i,j>0 



Definition 2.4. 



P 7 :=T( 7 )T( 7 - 1 ; 



y\(i+fc)( 7 ~i)(-i-fc) 



fc>l 



(4) 



»,j>o 



Thanks to the fact that 



EEii^ii^E^ii^ 00 



the product we have written on the right of (jlj) is a product of two Hilbert-Schmidt operators. 
So Vj differs from the identity by a nuclear operator. Hence its determinant is well defined 
(see for instance [28 J. In our notation we obtained the equality 



£00(7) = det(7> 7 ) 



(5) 



We will call "P 7 Plemelj's operator as it is related in a clear way with a Riemann-Hilbert 
factorization problem (see section 4) already considered by Josip Plemelj in 1964 [13j . 
Unfortunately, in concrete cases, det('P 7 ) turns out to be really hard to compute; nevertheless 
we can use some shortcuts also provided by Widom in his works ([6], [7] and [8]). 
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Proposition 2.5 (|_6j). Suppose that 7 satisfies conditions imposed in Szego-Widom theorem 

and, moreover, 7W = for i > j + 1 or 7W = /or i < j + 1. 

T/ien 

^00(7) = D^Miy (6) 

Proposition 2.6 (|8j). Suppose we have a symbol 7 satisfying conditions imposed in Szego- 
Widom theorem. Suppose moreover that 7 depends on a parameter x in such a way that the 
function x — ► 7(2;) is differentiable. If 7 -1 admits two Riemann-Hilbert factorizations 



such that 



1 - 1 (z) = t + {z)t_{z)=s-{z)s + {z) 

E 



fc>0 



fc>0 



«-(*) :=E 



fc<0 



fc<0 



Then 



trace 



({d x t + )t- - (5 z s_)s + ^9 a; 7 



(7) 



Also Dn(~i) can be expressed as a Fredholm determinant as pointed out for the scalar 
case in [9] and generalized for matrix case in [10J. 



Theorem 2.7 (Borodin-Okounkov formula, [ID])- Suppose that our symbol 7(2) satisfying 
conditions of Szego- Widom 's theorem admits two Riemann-Hilbert factorizations 



7(z) = 7 + (z)7_(*) = 0_(z)M*) 



such that 



7+ («) := E 7 - 
fc>0 

7 _(*) :=X)7^ 

fc<0 



fc>0 

m*) -E 



fc<0 



and G(7) = 1. Then for every N 

D N ( 1 )=D 00 (i)det(I-K 1 , N ) 

where, in coordinates, we have 

if min{i, j} < N 



(8) 



,Efcli(7-^ 1 ) (i+fc) (^- 1 7+) ( ~ j " fe) otherwise. 



Remark 2.8. One can easily verify that #I 1 7+ is the inverse of 7-0+ 1 so that, again, we deal 
with operators of type T(4>)T((f)^ 1 ) with = . Also we want to point out that the 

assumption (7(7) = 1 is not necessary. The formula for G(j) 7^ 1 is written in since in 
our case we will always have £7(7) = 1 we wrote the formula as it was given in |10j . 
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3 r functions for constrained KP and G elf and- Dickey 
hierarchies as block Toeplitz determinants. 

It is well known that given a point W £ Gr its time evolution 

W{t) := exp(tiz + t 2 z 2 + t 3 z 3 + . . )W 

is nothing but KP flow. Moreover points W £ Gr^ correspond to solutions of Gelfand-Dickey 
hierarchies (this is the celebrated Grassmannian formulation of KP hierarchy due to M.Sato, 
see for instance [3]). In their paper pQ Segal and Wilson gave a formula for the corresponding 
t function Tyy as determinant (Fredholm determinant) of the projection of W(t) onto H+. 
Here we take a slightly different approach that generalizes what has been done by Itzykson 
and Zuber in the study of Witten-Kontsevich r function in [15] (see also [TB] and [18] ) . This 
approach allows us to define not just T\y but also a sequence of {tw,n}n>o approximating 
rjy and being themselves t functions for some reductions of KP. 

Suppose we have an element W £ Gr^ n ^ ; thanks to results stated in Section 1 we can represent 
this element as 



W 



\Wi n . . 

with W{z) G LiGl(n,C). 

Also we assume that the matrix W(z) 



w nl \ 



■ ■ w nn J 



fff ) = W(z)H 



(n) 



W) 



{Wij{z)}iJ-_ 

--i + o(h) 



l..n 



satisfies 



Wij = z{0{\)),i > j 



w 



it 



, . (n) 

This means that we restrict to the big cell, i.e. we assume W = H\ . In the sequel we will 
always assume that W belongs to the big cell. Obviously we have a base for W £ Gr^ given 
by 

{z s Wj : s e N, j = l...n} 

where Wj is the column vector (wij-.-Wnj) 1 - . Using the isomorphim S : — > H the 
corresponding base for W £ Gr is given by 

{ujns+j = z ns Z{wj) : s e N, j = 1 . . . n} 

and, as in Section 1, we have 



i=l 



Note that thanks to the big-cell assumption we have u! ns+ j(z) = z ns+ i 1 (1 + O(i)). 

For these points W £ Gr^ n -* and vectors spanning them we define the standard time evolution 
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(KP flow) given by 



uj ns+j (t; z) : = exp 



. i>0 / 



Now we want to define the r function associated to W as limit for N — > oo of some block 
Toeplitz determinants twn- 



Definition 3.1. Take M = Nn a multiple of n. 



i~w,N(t) ■= det 



z l Uj{t;z)dz 



(9) 



l<i,j<M=Nn 



Fist of all we want to prove that tw,n is a block Toeplitz determinant and write explicitely 
the symbol. 

Lemma 3.2. For every j = 1 . . . n we have 

Wj(t;z) := E _1 (u;j(i,2:)) = exp(£(i, A))v)j(z) 

where we denote 



A 



/0 ... 

1 
1 



z \ 
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Proof We simply verify that multiplication by z on Gr corresponds to multiplication by A 
on Gr^ n ) through the isomorphism H — 1 . □ 

Proposition 3.3. twn is the N -truncated (n x n) — block Toeplitz determinant with symbol 

W(t;z):=exp(Z(t,A))W(z) (10) 

Proof Take i,j < n and s,v < N; the (i + sn,j + im)-entry of the matrix in the right hand 
side of Q is given by 



z- l - sn to j+vn (t- z)dz = j> z~ l ~ sn z vn uj j (t;z)dz = 

- i+{v - s) n Wjlit^Z^+^dz = Wji(t)^ 



,1=1. .n 

so that the right hand side of ([9]) is the transposed of the iV-truncated n x n block Toeplitz 
matrix with symbol W(i; z). □ 

In the sequel of this paper we will call such symbols Gelfand-Dickey (GD) symbols. 
Now generalizing what has been done by Itzykson and Zuber in [15] we expand Tw,N(t) in 
characters. In this way, assigning degree m to t m , it's possible to state a certain property of 
stability for tvi/,tv(£); namely every term of degree less or equal to Q will be independent of 
N for Q < N. 
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Proposition 3.4. 



rwMt) = e (n^i^ 



h,...,l nN >0 i 

where X = diag(xi, ...,x n N) is related to times {ti} through Miwa's parametrization 

(x k \ 

t k : = trace I — I 

and 



det 



/ x h+nN-l 


U+nN-2 

x x 


. . x x 


\ Ji+nN-1 


U+nN-2 


J-nN 


(xf- 1 


„nN-2 
x l 




det ... 




;■) 


\ r nJV-l 

\ x nN 


x nN 





Xh,...,l nN ( X ) 



Proof We start from determinant representation ([9]). The (z,j')-entry of the matrix will be 

2 w i i^h-iW 



where for every n > 



Pn(t) :=^-l ^ mt : Z)) dz 



2iri 



are the classical Schur polynomials and p n (t) = for every negative n. Then resumming 
everything we obtain 



TW,N(t)= ^2 (il^j J det b%+i-l(*)]ij=l. 
ki,...,k nN j 



.nN 



with /cj > 1 — j- 
Equivalently we write 



Tw,iv(i) = ^2 [Y[ u j lj+J det\pi j - j+i (t)]i ! j = i... nN . 
h,...,l nN >0 j 

On the other hand it's well known that under Miwa's parametrization this last determinant 
can be written as Xh,...,i nN (X) (see for instance [E],[16j); this completes the proof. □ 

Now it's easy to see that assigning degree 1 to every xi (which is equivalent to assigning 
degree m to t m ) we obtain 

M 



deg(x h ...i nN ) =J2 li 



i=l 



From this easily verified property we obtain the following 



11 



Lemma 3.5. Suppose deg(xh,...,l nN ) = Q < n N ■ Then, if the character is different from 
zero, we have 

Xl\...l nN XIi,:;Iq,0,...,0- 

Proof Suppose ^ 0, j > Q and k = V i > j. 

The j th column of the matrix [pi 3 -j+i(t)} has positive subscripts l\ + j — 1, 12 + j — 2, . . . , lj. 
On the other hand li = Q; hence the sum of these subscripts is 

i-i i 

r=0 r=0 

hence two subscripts must be equal, then two lines of the matrix are equal. □ 
From this corollary it follows directly the following result. 

Proposition 3.6. Up to degree Q the function T\y,N(t) does not depend on N with N > Q. 
This proposition allows us to define in a rigorous sense 

T W (t) := lim r WiN (t) (11) 

N— »oo 

where the limit is defined as limit of formal graded series in tf, this means that 

lim deg(T W:N (t) - T W {t)) = 00 

N— ♦00 

We will say that Twit) is the stable limit of Ty/,Nif) for N — > 00. 

On the other hand, in the sequel, we will prove that the symbol W(t) satisfies Szego-Widom's 
condition for every values of tj so that the limit in (jlip exist pointwise in time parameters 
and can be written as a Fredholm determinant. 

Now, following again [15] ; we write a differential operator A^jv(i) associated to the function 
Tw,N(t)- I n the sequel we will always write D for the partial derivative with respect to t±. 
We will prove that for every N the pseudo-differential operator Aw t ]\f(t)D~ nN satisfies Sato's 
equations for the dressing (see [H]) and we recover the usual relation between r and wave 
functions. 

Lemma 3.7. Define 

fs,N{t) : = J2^ k) Pk+nN-l(t) 

k>s 

Then we have 

rwAt) = Wr(/i,7v(t), • • • , fnN,N(t)) := det[D nN -> f iiN (t)]l<iJ<nN 
Proof From definition 13.11 the (i,j)-entry of matrix defining rw,N(t) is 

^U^pk+i-tit) 

k>j 

On the other hand we have 

D nN - j fi, N (t) = D nN -> (E-rWniV-i(t)) = ^^V'-iW 

V k>i ) k>i 

(using the equation D s {p m [t)) = p m - s (t)). Hence we obtained the proof. □ 
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Definition 3.8. We define the differential operator A\y,N(t) of order N in D as 

» rn _ Wr(f,f hN (t),...,f nN>N (t)) 
WMI} - Wr(f 1>N (t),...,f nN , N (t)) 

where f E depends in a differentiate way on {£j}i>i. 

Proposition 3.9. The following equations for time-derivatves of Aw,N(t) holds: 

JL(AwMt)) = ( A wM t )° iA WN( t ))+ " &w,N(t) Di (12) 

Proof It is enough to prove the equality of the two differential operators when acting on 
fi,N(t), ■■■fnN,N(t) which are nN indipendent solutions of the equation 

(A WjA r(i))(/(i)) = 

But this amounts to proving 



d 

— {A w>N ){t) 
9u 



f jtN (t) + A WjN (t)-^f j>N (t) = Vj 



which is true iff 

— (A WjN (t)f j , N (t))=0 Vj. 

This equality is obviously satisfied. □ 

Multiplying Apyjv from the right with D~ nN we found a pseudodifferential operator that, in 
fact, gives a solution of KP equations. 

Definition 3.10. 

Sw,N(t) ■= Aw,N(t)D~ nN 

Proposition 3.11. Sw,N{t) * s a monic pseudo-differential operator of order satisying Sato's 
equation (see for instance 

^-(SwMt)) = (SwMt)D i Sw 1 N ( t ))+ - SwMt)D l (13) 
Hence the monic pseudo- differential operator of order 1 

Cw,N(t) ■= (Sw,N(t)DS^ N (t)) (14) 
satisfies the usual Lax system for KP 

= [(£w,n)+'£w,n) (15) 

Otk 

Proof It is obvious that S\y,N is a monic pseudo-differential operator of order since A\y,N, 
which is of order nN, is normalized so that the leading term is equal to 1. Equation (|13p 
follows directly from (I12p . The derivation of Lax system from Sato's equations is well known: 
one has just to derive the relation 

£w,nSw,n = Sw,nD 
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for tf. and use the obvious relation [£w,N, £ k 



W,N\ 







□ 



It remains to prove that twn{^) is really the t function for these solutions CwN{t) of KP 
equations. We recall the usual relations between the dressing S, the wave function tp and r 
function given by 



V>(£; z) = S(i)(exp(£(i, z))) = exp(£(t, z)) 



r(t) 



(we recall that the notation t — A stands for the vector with i th component equal to ij 
All we have to prove is the following 



Proposition 3.12. 



ipW,N(t;z) := S W)N (t)(ex.p(£(t,z))) = exp(£(t,z)) 



T~W,N(t 



TW,N{t) 



(16) 



Proof Equivalently we prove that 



nN 



T~W,N(t 



T"w,Ar(i) 



Since we have 



Pn(t 



l) =Pn(t) ~ Z 1 p n -l(t) 



the right hand side of the equality above can be written as 

det 



D nN-l fi _ z -l D nN fi _ _ _ ^ _ ,-1^ 



z nN e £(x,t)_ 



JV-l 



/nTV — ^ 1 D nN f n w 



fnN - Z l Df n N, 



Wr (A,..., / nJV ) 

(here derivative is with respect to £1, we don't write dependence on /j on t to avoid heavy 



can 


be written as 




det 


/ z nN e ((t,z) 


• ••• /i 




\D nN f nN .. 


• • • • /niV / 




Wr (/!,... 


i /niv) 



It is easy to check that these two expressions are equal. □ 

We want now to study the structure of Cw,N with more attention; our investigation will lead 
us to discover that, actually, we are dealing with rational reductions (|19|.|18j) of KP. 
First of all we recall a useful lemma (proof can be found for instance in |29|). 

Lemma 3.13. Let {gi, g m } be a basis of linearly independent solutions of a differential 
operator K of order m. Then one can factorize K as 



with 



K = {D + T m )(D + T m _i) • • • (D + Ti) 

= Wr(gi, . . . ,g m -i) 
3 Wi(gu . . . , g m ) 
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We will state now properties of symmetry for / S) jv that will be useful in the sequel. 
Proposition 3.14. The following equalities hold: 

fs,N+l(t) = fs-nMt) (17) 

f s , N (t) = D n f s+n , N (t) (18) 

Proof We will use the equality 

which follows from the very definition of these coefficients. 
Then for (|17f) we have 

fs,N+l(t) = ^2^ k) Pk+n+nN-l{t) = ^U^n n) Pk+n+nN-l{t) = 
k k 

^2^s-nPk+nN~l(t) = f s -n,N(t) 
k 

For (fTBl) we have 

) s+n Pk+nN-l{t) - 2_ 
k k 

Y J ^~ n) Pk-n+nN-lit) = f SjN (t) 



D n (f s+n:N (t)) = D n ^2uj ( s+ k Jp k+nN „ 1 (t) = ^2uj ( s+ k Jp k+nN _ 1 _ n (t) 



□ 

Theorem 3.15. For every N > the pseudodifferential operator J~-w,N and its n th power 



can be factorized as 

• Cw,n = £i,w,n{£-i, w,n)~ 1 

• Lw,N = Mi t w,N(M 2i w,N)~ 1 

where all the factors are differential operators and 

• ord(Ci t w,N) = nN + 1, ord(£2,w,Af) = nN 

• ord(Mi ) w i jv) = 2n, ot&.{M2,w,n) = n 

Proof The first factorization comes directly from the fact that 

C W:N (t) = A T y iiV (t)D(A^iv(t))" 1 

For the second factorization we note that we have the factorization 

L w>N (t) = A WiN (t)D n (A W:N (t))- 1 

where the first operator A\v,N(t)D n has order M + n while the second (i.e. Aw,n) has order 
nN. Moreover as follows from (I18p we have 
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• A W)N f it N = Vi = 1, . . . , nN 

• A WjN D n f ijN = Vi = n + l,...,nN 

hence using lemma 13.131 one can simplify factorization above as 

Lw,N = Mi } w,n(M2,W,n)~ 
where M2,w,N is given explicitely by the formula 

M 2 ,w,N = (D + K n>N )(D + K n „ 1>N )...(D + K 1>N ) 

with 

/ Wr(/ n+ i,jy, 



K 



D 



log 



, fnN,N, fl,N, fj-l,N) 



\ Wr(/ n+1> jv, • • • , fnN,N, fl,N, ■ ■ ■ , fj,N) 



□ 



(19) 



Theorem 3.16. The sequence {£w,n}n>i satisfies recursion relation 

£w,N+i = TnC-WinCTn)' 1 
with 

T N = (D + T nyN ){D + r„_ 1)JV )...(D + T ljN ) 

rr m f^ T (fl,N,—,fnN,N,fl,N+l,—fj-l,N+l)\ 
Tj = D log — J - — . 

y WT{fi tN ,...,J n N,N,Jl,N+l,---Jj,N+l) J 
PROOF We observe that thanks to (fT7|) 

• A W)N f i>N = Vi = 1, . . . , nN 

• A W)N+1 f i)N = \/i = n + 1, . . . , nN 

• A W}N+ if itN+ i = Wi = 1, . . . , n 

Hence using again (|3.13f) we obtain the recursion relation 

Aw,n+i = T~nAw,n 

and from this last equation we recover the recursion relation for the Lax operator. □ 

Using the notation of [2D] we say that, for every N, tw,n £ cKPi in Ar n cKP„ in . 
This means that given a r function for KP with corresponding Lax pseudodifferential operator 
C we say that r G cKP m>n iff C m can be written as the ratio of two differential operators 
of order m + n and n respectively. These special reductions of KP begun to be studied in 
1995 by Dickey and Krichever (|18j.[19]); a geometric interpretation of corresponding points 
in the Grassmannian has been given in [21] and [22] . Going back to theorem 13.161 we point 
out that the first decomposition as well as the recursion formula are already known and, as 
pointed out in [20] . come simply from the fact that we have a truncated dressing. Actually 
our sequence of {tw,n}n>i is a part of a sequence already studied by Dickey in [17] under 
the name of stabilizing chain; in that article Dickey already provided the recursion formula 
written above as well as some differential equations for coefficients of T/v- Nevertheless, to our 
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best knowledge, connection with block Toeplitz determinants never appeared before. Also the 
fact that Tyy N £ cKP niI1 is something new. It could be interesting to find recursion relations 
as well as differential equations for M\ y w,N and Mz,w,N\ we plan to do it in a subsequent 
work. Till now all we can do is to infer from recursion formula for Lax operator the following 
formula 

(M\,w,N+l)~ TnM\ % W,N = (Mz^N+l)' ?nM2,W,N- (20) 

Now we want to go one step further and see what happens for iV — > oo. Obviously thanks 
to the property of stabilization stated in proposition 13.61 we can define a pseudodifferential 
operator Cw and a wave function ipw related to T\y in the same way as for finite iV and we 
will obtain a solution of KP as well. Actually a stronger statement holds. 

Proposition 3.17. Given W £ Gr^ the functions tw , ipw and L\y '■= {C-w) n ore respec- 
tively the t function, the wave function and the differential operator of order n corresponding 
to a solution of n th G elf and- Dickey hierarchy. 

Proof It is known pQ that subspaces satisfying z n W C W correspond to solutions of n th 
Gelfand-Dickey hierarchy. What we have to prove is that Lyy(t) = (Lw(t))+- 
In the sequel we will omit dependence on times of Ly/ and S\y if no ambiguity arises. 
From the usual relation 



we obtain immediately 

so that we have to prove that 



dt n 
dS w 



(t;z) = (L w (t))+^ w (t;z) (21) 
+ SwD n = {L w ) + Sw 



9t n 

dSw 

dt n 



On the other hand 



ipw(t; z) = exp(£(i, z)){\ + ^2 s i(*) 2 ^ 

1=1 



where 

oo 



S w (t) = l + ^ Si {t)D- 1 
1=1 

Using this explicit expression for the wave function and substituting in (|21j) we obtain 

dsi{t) 



(L w (t)) + ip w (t; z) - z n ^ w (t; z) = exp (£(t, z)) ^ ■ 



dt r . 



-z 



The left hand side of this equation lies on W(t) = exp(^(t, z))IU for every t so that multiplying 
both terms for exp(— z)) one obtains that they belong to subspaces transverse one to the 
other (W and H_), hence both of them vanish. This means that §^ = for every i. □ 

In virtue of this proposition, when computing tw associated to W G Gr^ n \ we will always 
omit times tj n multiple of n. Setting {tj n = 0, j G N} will be important in order to be able 
to apply Szego-Widom's theorem; in this case we will write t instead of t. 
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Proposition 3.18. Take any W G Gr (n) in the big cell of Gr* n > and a corresponding GD 

symbol W(t; z). 

Then 

Tw(t) = det(V w( t; Z ))- (22) 

Proof All we have to prove is that conditions of Szego-Widom's theorem are satisfied and 
G(W(t; z)) = 1. We observe that 

W(t;z) G LiGl(n,C) Vt 

2 

since we can always find W(z) G LiGl(n, C) such that W = W(z)H^ and exp(£(i, A)) is 

continuously differentiable (obviously when restricted to a finite number of times) . 
Moreover 

det[exp(£(£, A))] = 1 

since we deleted times multiple of n and det(W(z)) = 1 + 0(z^ 1 ) by big cell assumption. 
This implies that we have 

A det (Wit; e ie )) = 

and 

G(W(t;z)) = l. 

□ 

We are now in the position to state the main result of this paper, proof follows from results 
obtained above. 

Theorem 3.19. Given any point W G Gr^™^ and corresponding GD symbol W(t;z) the 
following facts hold true: 

• {rw,N(t) := D N (W(t; z))} < N<oc 

is a sequence of t functions for KP associated to wave function i(^w,N(t', z ) o,nd 
pseudodifferential operators Cw,N{~t) given respectively by ( fifil) and |i^| ), 

• For every N > we have tw,n G cKPi )n Ar H cKP„ jn . 

Explicit factorizations of Lax operator and its n th power are given in theorem \3.15l 

• The sequence admits stable limit Tw(t)- 

The stable limit is a solution ofn th G elf and- Dickey hierarchy. It is equal to the Fredholm 
determinant det Vyya. z >. ■ 

Remark 3.20. Also all the Tw,N(t) can be expressed as Fredholm determinants; in order to 
give explicit expressions we need a certain Riemann-Hilbert factorization of symbol W(t; z). 
This factorization will be obtained in section 5 and it will be exploited to express Tw,N(t) as 
a Fredholm determinant. 
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4 Riemann-Hilbert problem and Plemelj's integral formula. 



It is evident from proposition 12,61 that Riemann-Hilbert decompositions of symbol 7 for a 
block Toeplitz operator plays an important role in computing ^00(7)- 

Here we will show that actually Plemelj's operator itself enters in a integral equation (see 
|13j ) giving solutions of Riemann-Hilbert problem 

( p + (z)= 7 T (z)^(z). (23) 

Here ip+(z) and (f~(z) are respectively analytical functions defined inside and outside the 
circle. In this section we consider a smaller class of loops; 7(2) will be a matrix-valued 
function that extends analytically on a neighborhood of S 1 . For convenience of the reader we 
recall here main steps to arrive to Plemelj's integral formula [13]. 

Lemma 4.1. Suppose that f+(z), f—(z) are functions on S 1 satisfying 

i/(C2)-/(Ci)i<ic 2 - arc- 

for some positive constants fi, C and for every Ci, C2 £ S . Necessary and sufficient conditions 
for f+{z) and f-(z) to be boundary values of analytic functions regular inside or outside 
S 1 C C and with value c at infinity are respectively 

1 I MO - /+(£) d( = (24) 



2ni J ( 
1 fUC)-f-(z) 



2m J C 



d( + f-(z)-c = (25) 



We have to point out that here both £ and z lies on S 1 so that one has to be careful and 
define (I24p and (I25p as appropriate limits. Namely one proves that taking £ slightly inside 
or outside S 1 along the normal and making it approach to the circle we obtain the same 
result which will be, by definition, the value of our integral. Now suppose we want to find 
solutions of ()23|) : we normalize the problem requiring ip_ taking value C at infinity. Taking 
an appropriate linear combination of (|24p and ([25]) and using (f2"3"P we find that <f-(z) must 
satisfy the equation 

Note that here we do not have to take any limit since the integrand is well defined for 
every point of S . We also want to consider the associate homogeneous equation 



as well as its adjoint 



Obviously, as usual in Fredholm's theory, the equations (|27p and (|28p either have only trivial 
solution or they have the same number of linearly independent solutions. 



o = «,) + -L / tMx^KW ( ) (28) 

ZTTl J Q — Z 
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Lemma 4.2. Consider two adjoint RH problems 

<p + (z)= 1 (z) T V-(z) (29) 
1> + (z)=<y(z)il>-(*) (30) 

normalized as ip-(oo) = i/?_(oo) = 0. 

Any solution ip^ of 129\) is a solution of {2H\) as well as any solution ip + of $30\) is a solution 

ofm- 



Proof We just repeat computations made for non- homogeneous case. □ 

Now we introduce a new ii 
to the Plemelj's operator. 



(n) 

Now we introduce a new integrable operator acting on H\ and prove that it is actually equal 



(n) 

Definition 4.3. For every f £ H\ we define 



[V^](z) := pr + L(z) + 7(Z)7 ( 7^ ) - mdC, ) (31) 







where pr , denote the projection onto ii+ . 



Proposition 4.4. 



Proof We write P~ in coordinates and verify we obtain the same as in dH). To do so as in 
the definition of integrals (|24|) and (I25p we compute (I3ip imposing |£| < \z\; the formula will 
hold when £ approach to S* 1 in the same way as in (|24p and (|25p . For a consistency check we 
will prove we obtain the same result imposing \(\ > \z\. Let's start with |£| < \z\; we have 

+ /e5('- £ 7<"(7- 1 )<«^c)E* ( "< 'f 

7 fc>i p.gez s>o * 

Imposing /c + g + s = 0we get that this is equal to 

pGZfe>ls>0 teZ k>\ s>0 

Taking the projection on H^l' we obtain exactly formula (jlj). Now for |£| > \z\ we have 



*w + i / E S ( E ^ (p) (t- 1 )^^ - /) £ ^ (s) c 



/c>0 p,<?€Z s>0 



c 



20 



Imposing q + s = k we arrive to 

^ 7 (p)( 7 - i )( fc - s ) 7 /,w^+p = ^ ^ 7 (t-fc)( 7 -ij(fc--)^,w^ 

k,s>opez k,s>otez 

Taking the projection on we obtain that this is equal to T(7)T(7 _1 ) so that the two 
computations for \(\ < \z\ and for \(\ > \z\ coincide in virtue of lemma 12.31 □ 

Theorem 4.5. Suppose we are given a symbol 7(2) analytic in a neighborhood of S 1 and such 
that 

^oo(7)/0 

Then the Riemann-Hilbert problem 

normalized as <^_(oo) = C admits (if existing) a unique solution. 

Proof Suppose we have two distinct solutions (y>j_, ^1+) and (<p2-, ¥>2+); taking the dif- 
ference we obtain a non-trivial solution of (129j> . Then also (|30p admits non trivial solutions 
and the same holds for Q28|) . But this means that we have a non zero ^(2) 6 -H+ such that 
[Vjip] (z) = which is impossible since 

det(7> 7 ) = Doo(7) / 

□ 

Existence of factorization will be treated in the next section for the specific case of Gelfand- 
Dickey symbols. For a general treatment of the problem of existence see [13] . 



5 Factorization for Gelfand-Dickey symbols 

Here we will prove that for Gelfand-Dickey symbols we can write the unique solution of 
factorization ([2*3*1) in terms of data L\y(t), ipw(t; z). We recall that Lyy(t) and i/jw(t;z) are 
the stable limits of Lw,N(t) and tpw,N(t', z). They represent the differential operator and the 
wave function associated to the solution Twit)- Our exposition here is closely related to [13] . 
At the end of the section we will use the factorization obtained to express any tw,n$) as a 
Fredholm determinant. As we have written before in the proof of proposition 13.171 we have 
the relation 

L w (t)ij w (t; z) = z n 4) W (t; z) (32) 
where ipw(t> z ) admits asymptotic expansion 

ip w (t; z) = exp(£(t, z))(l + 0(z~ 1 )) 

Now out of ipw w e construct n time-dependent functions 

^W,S\ z) ■= D l (ip w (t; z)) : i = 0, . . . ,n- 1 

belonging to the subspace W G Gr. 



21 



Definition 5.1. 



/i Ci • 

1 C2 • 
\l Cn • 



n-1 



/^w,o(*;Ci) ^w,i(£;Ci) 
Vw,o(*;C2) V , w,i(*;C2) 



• • CT7 W,o(t;Cn) ^w,i(t;Cn) • 

where Q is the i th root of z. 

Proposition 5.2. The matrix ^w(t; z) admits asympotic expansion 

*w(i;A) =exp(e(t;A))(/ + 0(^ 1 )) 
Moreover under the isomorphism S _1 : H — > we can write W S Gr^"^ as 

W = * w (0,z)hW 



ip w ,n-i(t; Ci)\ 

l(lW,n-l(t;Cn)J 



(33) 



Proof One has to note that the i th column of matrix ty\y(t; z ) is nothing but H 1 {ipw,i(^ z )) 
so that asymptotic expansion follows easily. Equation (|33p corresponds to the fact that 
{z ns *p w ,i(0, z) : s G Z} is a basis for W. □ 

Observe that, since we also have 



W = W(z)H. 



we obtain 



*w(0,z) =W(z)(I + 0{z- 1 ). 
From this equation and from lemma [231 it follows that for every A r >0we have 

T N (w(t; z)(l + 0(z- 1 ))) = T N (W(t; z))T N {I + O^" 1 )). 

Now since for every iV 

detCTjvCJ + O^- 1 ))) = 1 
we will assume, without loss of generality, that 

* w (0,z) = W(z) 

since this is true modulo an irrelevant term that does not affect values of determinants we 
want to compute. We now want to define a matrix <&\y(t;z) analytic in z near and with 
similar properties as i&w(t;z). 

Definition 5.3. Let <f)w(t;z) be the unique solution of 

L w (t)<p w (t;z) = z n <p w (t;z) 

analytic in z = and such that 

(Z>V)(0,z) =z { :i = 0,...,n- 1 
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We define 



/i ci ••• cr\ 



*w(t;z) := 



i C2 ... C 



n-l 



/0w,o(*;Ci) 0w,i(*;Ci) 
4>w,o(t;(2) (f>w,i(t; (2) 



<j>w,n-i{t; 
0w,n-i(*; C2 



\1 Cn ••• C V \<£w;o(*;Cn) 4>w,i(t;Cn) ■■■ <j>w,n-i (i; Cn)J 
where as before Q is the i th root of z and 

(pwAt) : = Di (<t>w(t; z)) ■ i = 0, • • • , n - 1. 

Remark 5.4. $w{t] z ) admits regular expansion in z = and Cauchy initial values we imposed 
on (j>w imply 

$ w (0;z)=I. 

Proposition 5.5. ^w(t] z )&w$'> z ) does not depend on ti for any i. 
Proof It is well known that equations 

d ± 
— / = (L^) + f 

satisfied by 4>w and ipw can be translated into matrix equations 

d 



Ot, 



F = FM 



satisfied by ^w{h z ) and $w(t', z ) ( one ca n write explicitely M in terms of coefficients of 

i_ 

(Lyy)+)- Hence we have 

<$w(t; z)^(t; z)<Z> w (t; z)M^(t; z) = 

Theorem 5.6. Given a Gelfand-Dickey symbol 

W(t;z) =exp (£(i,A))wO) 

one can factorize it as 



□ 



W{t;z) 



exploit, A)W w (-t,z) 



where the term inside the square bracket is analytic around z = 00 and the other is analytic 
around z = 0. For assigned values of t for which 

T W (t) + 

this is the unique solution of the factorization problem normalized at infinity to the 
identity. 
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Proof Using the previous proposition we have 

W(t;z) = exp£(t, A)W(z) = exp(£(t,A))* w (p,z) = 
exp(f (t, A))*(-t; 2)*^(-*l ^)^w(0; *) = exp(£(t, A))*(-t; 2)*wH; z) 

Unicity of the factorization follows from section 4. □ 
Corollary 5.7. For every N > 

nv,jv(*) = det(7 - ifw(tj^),jv) 

with 

i/ min{i, j} < N 



W(t;z),N)ij 



X£i(*VK(-t;z))( i+fc )(^ H /(-t;z)- 1 )(^- fe ) otherwise. 



Proof It is enough to apply Borodin-Okounkov formula using factorization obtained above. 
□ 



Corollary 5.8. For every N > 

TW,N(t) 



det [T^wi-t^z^T^wi-t-z)- 1 



T W ,N+l{t) ^ ; Jn,N 

(observe that the right hand side of this equation is an ordinary n x n determinant, not a 
Fredholm determinant). 

Proof 

rw,N(t) _ det ( J - k w(p,z),n) 



TW,N+l(t) d^i 1 ~ K W(t;z),N+l) 

On the other hand the operator (I — Ky^a.^ jv + i) — K w a. z \ N ) can be written as a block 
matrix obtained taking the identity matrix and replacing the N th block column by the N th 
block column of the matrix with (i, j)-entry equal to 



j^^wi-t-z^^wi-^z)- 1 : 



1 vv \ v 1 vv \ 

k=l 

Hence proof is obtained applying lemma 12.31 □ 



6 Rank one stationary reductions 

and corresponding Gelfand-Dickey symbols 

We want to describe, more explicitely, GD symbols corresponding to solutions of Gelfand- 
Dickey hierarchies obtained by rank one stationary reductions. In order to emphasize that 
we are dealing with rank-one generic case instead of the standard expression Krichever locus 
we will speak about Burchnall-Chaundy locus. 
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Definition 6.1. Given a point W G GA n ^ we say that W stays in Burchnall-Chaundy locus 
iff the Lax operator L\y of the corresponding solution satisfies 

[L W ,M W ] = 

for some differential operator M\y of order m coprime with n. Without loss of generality we 
also assume m > n. 

The name we use is due to the fact that, already in 1923, Burchnall and Chaundy were 
the first to study algebras of commuting differential operators in [23] where they stated this 
important proposition we will use in the sequel. 

Proposition 6.2 (|23j). Given a pair of commuting differential operator L, M with relatively 
prime orders it exists an irreducible polynomial F(x, y) such that 

F(x,y)=x m + ...±y n 

and F(L, M) = 0. 

This proposition in particular allows us to associate to every Burchnall-Chaundy solu- 
tion a spectral curve defined by polynomial relation existing between the pair of commuting 
differential operators. From the Grassmannian point of view one can define an action A of 
pseudo differential operators in variable t\ on H by 

A : ^DO x H — > H 
and, using this action, prove the following propostition 

Proposition 6.3 ([!]). Given a point W in the Burchnall-Chaundy locus one has 

z n W C W (34) 
b(z) C W (35) 

where L\y and M\y are of order n and m respectively and b{z) is a series in z whose leading 
term is z m . Conversely, if W satisfies above properties, it stays in the Burchnall-Chaundy 
locus. 

Proof We just sketch the proof and make reference to Mulase's article Suppose we are 
given Lyy and My/'-, under conjugation with the dressing Swif) we have 

S^(i)L w (t)S w (t) = 

Under the action A this gives invariance of W with respect to z n while invariance with respect 
to b(z) is obtained acting with 

S^-(t)M w (t)S w (t) 

Viceversa given W we reconstruct the dressing Sw(t)', using it we define Ly/(t) and My/{i) 
conjugating pseudodifferential operators corresponding to z n and b{z). In particular observe 
that also z n and b{z) will satisfy the same polynomial relation as Lw(t) and My/(t). □ 
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Remark 6.4. Without loss of generality we can assume 

hj^ dz={] VseZ - (36) 

Now suppose we are given an element W = W(z)H^ 6 Gr^ in the Burchnall-Chaundy 
locus. Using the explicit isomorphism H we can construct a matrix B(z) := 6(A) such that 

B(z)W C W. (37) 

Proposition 6.5. 

C(z) := W- 1 (z)B(z)W(z) 

has the following properties: 

• C(z) is polynomial in z. 

• trace(C(z)) = 

• m = maxj(j — i + n deg Cij(z)) Vj = 1 . . . n 

• T/ie characteristic polynomial Pc(z)W °fC(z) defines the spectral curve of the solution. 
Proof Equation (|37p can be equivalently written as 

W- 1 {z)B{z)W{z)H { ^ ) C H { + ] 

and this means precisely that C{z) can't have terms in z~ k for any A; > 0. 
The other properties are satisfied if and only if they are equally satisfied by B(z) so that 
we will prove them for B{z) instead of C(z). B(z) is traceless thanks to equation (|36|) and 
thanks to the fact that 

trace (A fc ) = \/k ^ sn 

The third properties is satisfied as B{z) = 6(A) represents in H multiplication by a series 
whose leading term is equal to m. For the last property we observe that if F(x, y) is the 
polynomial defining the spectral curve, i.e. F(Lyy, Mw) = 0, then we will have 

F(dwg(z,z,...,z),B(z)) = 

as well; on the other hand thanks to Cayley-Hamilton theorem we have 

PB(z)(B(z)) = 0. 

Since F is irreducible and Pb( z ) W nas the same form 

^(,)(A) = A" + ...±z m 

we conclude that they are equal. □ 

Observe that since W(z) is defined modulo multiplication on the left by invertible trian- 
gular matrices also C{z) is defined modulo conjugation by elements of the group A of upper 
trianguar invertible matrices. It was a remarkable observation of Schwarz [23] that actually 
Burchnall-Chaundy locus can be described by means of matrices with properties as in propo- 
sition E3] modulo the action of A. Here we adapt the results of [21] to our situation. Namely 
we explain how, given C(z), one can recover W(z) and the corresponding spectral curve. 
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Proposition 6.6. Given a matrix C(z) such that: 

• C(z) is polynomial in z. 

• trace(C(z)) = 

• m = maxj(j — i + n deg Cij(z)) Vj = 1 . . . n 

it exists a unique W = W(z)H^ in Burchnall-Chaundy locus such that its spectral curve is 
defined byp C ( z ){\). 

In order to prove this proposition we need two lemmas. 

Lemma 6.7. Given a polynomial matrix C{z) such that 

m = max(j — i + n deg Cij(z)) Vj = 1 . . . n 

i 

(with m and n coprime) coefficients of characteristic polynomial 

Pcw(A) := A™ + ci(^)A n - 1 + . . . + c n (z) 



satisfy 



n deg c s < ms Vs = 1, . . . , n — 1 
deg c n = m 



Proof From 



n deg Cjj < m — j + i 
and definition of determinant follows immediately that 

n deg c s < ms Vs = 1, . . . , n. 

Strict inequality for s < n follows from the fact that m and n are coprime. For the equality 

deg c n = deg^det(C(z))^ = m 

we observe that in every line there is a unique element Cij(z) such that m = j — i+ndegCij(z); 
taking this unique element for every line and multiplying them we will obtain the leading term 
of determinant which will be of order m. □ 

Lemma 6.8. The equation 

\ n + ci (z)X n ' 1 + ... + Cn(z) = (38) 

with 

n deg c s < ms Vs = 1, . . . , n — 1 
deg c n = m 

and n, m coprime has n distinct solutions {Aj = b(Q), i = 1 . . . n} with 

b(() = C(i + o(c 1 )) 

(as usual Ci is the i th root of z). 
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Proof Imposing Aj = £|™ we have a solution of the equation 

(C) n + ^CXCrT" 1 + • • • + c n (C) = o 

at the leading order ran. Then imposing Aj = CJ"^ + ^iC rl ) an d plugging it into the equation 
(1581 one obtains 



(cr + JicrT + ci(cd(c + hcr~T~ L +■■■+ = o(cr n ) 

Zi can be found so that terms of order nm — 1 in the equation vanish; going on solving the 
equation term by term we obtain 

Clearly coefficients lj do not depend on the choice of the root Q so that it exists 6(A) with 
stated properties. □ 



Now we can prove proposition 16.61 

Proof We start computing the characteristic polynomial Pc(z){^)\ thanks to lemmas 16.71 and 
16.81 we find n distinct roots 6(Ci ) , • • • , b(Q n ) with properties stated above. 
The aim is to find W(z) such that 

W(z)b(A)W- 1 (z) = C{z) 

Since we have n distinct solutions {b(Q), i = 1, . . . , n} of the equation 

PC(z)W = 
it exists a matrix T(£i, . . . , Cn) such that 



T(Ci)C(z)T" 1 (C<' 



6(Ca) 


V o ... 



o \ 




b((n)J 



On the other hand it's easy to observe that the matrix A can be diagonalized as 



A 



and this means that multiplication by b(z) can be written in as multiplication by 



/I 


Ci • 


■ cr'\ 




fCi o . 


.. 0\ 




( i 


Ci • 


• • cr x \ 


1 


c 2 • 


f-n—l 
■ ■ ^2 




o C 2 • 


.. 




l 


C2 • 


fTl— 1 










o ... ' 


•• 












Cn ■ 


•• Q- 1 ) 




\o ... . 


•• Cn) 




v'i' 


Cn • 


• • crv 



/i 


Ci • 


• • crv 


- 1 /6(Ci) o . 


• \ 


/I 


Ci • 


•• cr x \ 


i 


c 2 . 


An— 1 
• • S2 


6(Ca) • 





1 


C2 • 


/•n— 1 
• • ^2 








... ' 













Cn • 


• • cr 1 / 


V ... . 


• KCn)J 




Cn • 


• • cr 1 / 
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Hence we have 





(i 


Ci • 


■ ■ cr 




w(z) = r- 1 (( l ) 


1 


C2 . 




-i 






Cn 


■■ Q~ 


V 



Note that W(z) is defined modulo the action of A so that, by construction, C(z) corresponds 
to a unique W G Gr^ n ^ such that 

W = W{z)H i ^ ) . 

□ 

Remark 6.9. As it was pointed out by Schwarz [23], matrices C(z) with properties stated 
above can be used to describe points in the Grassmannian describing string solutions of 
Gelfand-Dickey hierarchies, i.e. solutions associated to reduction of type 

[L, M] = 1 

This class of solutions has not been treated in this article since they do not live in Segal- Wilson 
Grassmannian but just on Sato's Grassmannian constructed on the space of formal series; this 
means that we cannot use any more Szego-Widom theorem as the analytical requirements are 
not satisfied. Nevertheless some results obtained in section 3 still hold since the property of 
stability for {rw,N(t)} does not depend on analytical properties of the symbol W(z). Hence 
one can try to apply the approach used in this article to the study of these (much less studied) 
string solutions; perhaps results obtained by Okounkov and Borodin in [9j for formal series 
and a generalization to block case can play in the setting of formal theory the same role played 
by Szego-Widom theorem in this paper. 

Example 6.10 (Symmetric n-coverings) . Take a symmetric n-covering C of P 1 given by equa- 
tion 

nk+l 

A" = P(z) = J! (z - a,) (39) 
j'=i 

For this particular type of curves, choosing in a appropriate way the divisor on the curve, we 
can write explicitely W(z),B(z) and C(z). We start to observe that for any W corresponding 
to this spectral curve we have b{z)W C W with 

b(z) =P(z n )^ 

Then it's easy to prove that the corresponding B(z) = 6(A) can be written as 

/ n-l , , 1 \ 

/ ... z— P(z)n\ 

Z~nP( Z )n ... 



\ ... z~P{z)n / 
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Now we define n functions 



i = 1, 



We take 

W := diag(u>i(z), ...,w n (z)) 
It is easy to verify that the matrix 

C{z) = W- l {z)B{z)W{z) = 



,n. 



( 












v ' w\(z) 



\ 







z ~ip{z)\ : ^4*l 

V / Wn(Z) 







/ 



is polynomial in z. It is worth noticing that this example already gives all possible double 
coverings; hence for any (possibly singular) hyperelliptc surface we found (assigning a par- 
ticular divisor) the GD symbol of the corresponding algebro-geometric rank one solution of 
KdV. 

Example 6.11 (Rational solutions). As pointed out by Segal and Wilson [lj, subspace of 
Burchnall-Chaundy locus corresponding to rational curves are given by W = W(z)H { + ] with 
W(z) rational in z. In particular the corresponding Gelfand-Dickey symbol will satisfy hy- 
pothesis given in proposition 12.51 so that we recover the following (known) result. 

Proposition 6.12. Every rational solution of Gelfand-Dickey hierarchies can be written as a 
finite-size determinant. 

For instance, for n = 2, taking 

'1-cPz- 1 



W{z) 



1 - c 2 z~ l 



HI 



the inverse of Gelfand-Dickey symbol is equal to 

W-\t;z) 



( cosh(zU^2 i>0 t 2 i + iz 2i ] 



-Z2 



2 sinh(z2 ( £ t2i+iz 2i ))\ 



z 



\-z 2 S inh(z2^.> t 2i+1 z 2i )J cosh(z2^.> t 2m z 2i )J j 



z-d? 




\ 



z — c 2 / 

Simply taking the residue one obtains that the corresponding r function will be equal to 

( coshfc^+i^ 1 ) -dsmh(Zi> t2i+id 2i+1 )\ 



Tw(h,t 3 , • ••) = det 



\-c sinh ^ X^i>o t 2 i+ic 2l+1 ) cosh ( £\> t 2i+ ic 
and recover 2-solitons solution for KdV. 



2i+l 
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We want to point out that, for algebro geometric solutions treated in this section, the 
problem of factorization for Gelfand Dickey symbol can be easily translated into a Riemann- 
Hilbert problem on some cuts on the plane with constant jumps. For simplicity we reduce 
to the case n = 2; the procedure used here is equivalent to the one used by Its, Jin and 
Korepin in [25] and generalized by Its, Mezzadri and Mo in [26]. Suppose we want to solve 
the factorization problem 

W(t; z) := exp fat, A))w(z) = T_(i; z)T + {t; z) 



for our GD symbol with W(z) = diag(wi(z), W2(z)) as in example I6,10| since it will appear 
many times we denote A the matrix 



A 

Also we impose 



i 

1 -Jz 



2g+l 

P(z) := H (z - a,) 
with all aj having modulo less then 1 and 

IKH < IKH < • • • < ||02 ff+ i|| 

We denote li,...l g+ \ the oriented intervals (ai, 02), (03, 04), ...(03^+1, 00). Instead of looking 
for T_(t; z) and T + (t; z) we define a new matrix S(t; z) imposing 

' S(t; z) := A exp ( - 0, A)) T_ (t; z) z>\ 

< 

k 5(t; z) := iW(^)T; 1 (?; z) z<\ 

Proposition 6.13. S(t;z) has the following properties: 

• It has no jumps on S 1 

• It has jumps on intervals lj; precisely calling SL(t;z) and Sji(t;z) the values of S(t;z) 
approaching from the left and approching from the right the interval we have 

S L (t;z) := K S R (t;z) 

• It is invertible in any points but aj; there it has singular behaviour of type 

'1 



with Sj(t; z) invertible in aj; minus is for a\, . . . , a g , plus for the others. 
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• At infinity it behaves as 
S(t;z)r 



exp y - y/z(tiz + t 3 z 4 
exp i^fz{t\z + t 3 z + 



y/zexp ^ - yfz{t\Z + t 3 Z+ . . .] 

—y/zexp ( yfz(t\Z + t 3 z + . . .) 



Proof Let's call £+(£; z) and SL (t; 2) the limiting values of S(t; z) approaching the unit circle 
from inside and outside; we have 

S + (t;z)SZ 1 (t;z) = Aexp{-^(t,z))T^(t;z)T + (t;z)W- 1 (z)A- 1 = 
Aexp(-^(t,z))exp(^,z))W(z)W- 1 (z)A- 1 = I 

and this proves we haven't any jumps on S 1 . 
Writing explicitely S(t;z) as 

/ exp ( - y/z{t\z + t 3 z + . . .) ) y/zexp ( - yfz{t\z + t 3 z + . . .) )\ 



Sii-z) 



S{t-z) 



^ exp (^/z(tiz + t 3 z + . . .)J -y/zexp + t 3 z + . . .)j j 



( 1 {P{z))^ \ 

nf= ( z - a o) 



T_(t;z) z>l 



Tl^fcz) z<l 



V n?=o(*-w 



we obtain almost immediately the other points of the proposition; the only thing we have to 
observe is that both T + (t;z) and T_(t;z) are invertible inside and outside the circle respec- 
tively. This is because we have 



det W(t; z) 



This condition combined with 



gives 



det(7+(t; z)) det(T_(t; z)) 



lim det(T_(t;z)) = 1 

2^00 

det(T+(i;z)) = 1 

{P{z))\ 



det(T_(t;js)) 



nLo( z - °i 



□ 



The Riemann-Hilbert problem given by propostion 16.131 is equivalent to the one proposed 
in section 5. What can be done is to write explicitely the solution S(t; z) using 6 functions 
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associated to the curve; this is what has been done in |25] and [26]. Actually comparing 
previous proposition with results obtained in section 5 we immediately realize that 



S(t;z) 



so that all we have to do in our case is to write down Baker- Akhiezer function in terms of 
special functions. We can carry on the same procedure for n arbitary; the only difference 
will be that the jump matrices will remain constant but more complicated; in any case the 
solution of this Riemann-Hilbert problem with constant jumps will be 



S(t;z) 



Vw,n-i(£;C2) 



WW,o(£;Cn) 4>W,l(t;(n) ■■■ i>W,n-l(t;(n)J 



Explicit formulas involving special functions can be used here to apply proposition 12.61 to our 
case. For instance taking the elliptic curve C given by equation 



w 2 = 4z 3 - g 2 z - g 3 



with uniformization given by the Weierstass p function 

(z,w) = (p(u),p'(u)) 

one can write wave function as 

a(u — c — x)a(c 

ip(x,t,u) :=- 



exp (x((u) - -tp'{u) 
a(u — c)a(x + c) \ 2 

(here £ and a are Weierstrass £ and a function respectively, x and t correspond to the first and 
the third time). With some tedious computations, making the change of variables u = u(z), 
the right hand side of equation ((Jj) can be obtained. It turns out that the only relevant 
factorization is the one given by 

W~ 1 (x,t;u) = W~ 1 (u)^(-x,-t;n) ^(-x, -t; u) exp(-xA - tA 3 ) 

where as before (we just wrote function of u) we have 

/ i_ \ — i 

/ 1 I (r\ ( 1 1 I I 0. \ 

^(x,t;u) : 



(p(u))3 
"(P(«)) 



t/j(x,t,u) d x i[)(x,t,u) 
ip(x,t,—u) d x ip(x,t,—u) 



Plugging into equation ([7]) we obtain 



d 



— r(x,t) = Kt + 2C(-c) - 2((x - c) 
ax 

(here K is some constant); taking another derivative we obtain elliptic solution of KdV as 
expected. 
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